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Abstract 

In this letter we explore whether we have a chance to observe a flavor-changing effect in r 
appearance experiments ICARUS and OPERA. 
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The atmospheric neutrino anomaly |l] and the solar neutrino deficit are well described 
by the neutrino oscillation. Therefore many experiments are being planed to observe the 
oscillation directly and some of them are now carried out |J. For example, the atmospheric 
neutrino anomaly is very naturally explained by the — > v T oscillation and for the direct 
measurement of this oscillation (transition) there are two experiments proposed at CERN, 
ICARUS and OPERA || , that aim at observation of tau neutrino appearance events and 
they will start to take data in 2005 f| [7j]. 

Moreover it is expected that almost all the oscillation parameters will be determined 
by oscillation experiments in near future (see for example ||). Though in studies of those 
experiments their main concerns are how well they can determine the oscillation parameters, 
it is also pointed out that a long baseline experiment can probe new types of flavor changing 
interactions with the oscillation manner 10, 11, [TJ]. Among such new effects to probe 
the fi-T flavor- changing effects is most important since in many models to explain a large 
mixing for <-> v T oscillation, large /x-r flavor changing interactions are accompanied. In the 



previous we investigated the feasibility to observe such exotic interactions in oscillation 
experiments and showed that the — > v T channel works most effectively to explore the 
yU-r flavor violating interactions |16|]. Fortunately, ICARUS and OPERA use this channel 
and hence in this work we examine the performance of ICARUS and OPERA for the new 
interaction search in the three-generation framework. 

First we briefly review the key idea for exotic interaction search in an oscillation experi- 
ment. In the following, we consider only the /i-r flavor violating effects. In a long baseline 
experiment, what we really observe are the signals caused by the secondary charged par- 
ticles such as muons. That is, we do not observe neutrinos themselves. In other words 
neutrinos are unobserved intermediate states. Therefore, if there are some kinds of new 
interactions that can induce the completely same final states as the standard model does, 
then the interference between these two amplitudes takes place. It means that the effect 
of new interactions appears with the strength of not the square of the exotic coupling but 
itself. 

The contribution of the new interactions can be divided into three stages, neutrino beam 
production, its propagation, and its detection. First we refer to the production process. 
The neutrino beam of the CNGS M facility is produced by pion decay. In this case we 
can parametrize the effect of the exotic interactions in pion decay as the shift in the muon 
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neutrino state from pure flavor eigenstate, \u a ), to mixed state, \u^) ||, because the helicity 
states of the muon and the neutrino are fixed in this decay: 

K) = U' m \u a ), f^ = (0,l,e* T ) (1) 

where e^ T denotes the ratio of the coupling of the exotic decay of pions, n + — > (i + v T , to the 
standard one. Next, during the beam propagation from CERN to GranSasso, the neutrinos 
feel the matter effect due to not only an ordinary interaction but also the new ones, which 
can be interpreted as the shift of the potential fpTOfl : 
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where a = 2y2GFn e E v is the standard matter effect whose origin is the weak interaction, 
n e is the electron number density, and E v is the neutrino energy. Finally at the detection 
process neutrinos are also affected by the new flavor-changing interactions. To parametrize 
these effect, in principle, will be very complicated since we have to deal with the hadronic 
process. However, in this letter, we assume for simplicity that we regard the neutrino state 
at the detection process as the flavor mixed state, Q, just like source state, 



K> = t£>«>, ^™ = (o,< r ,i). (3) 

If more precise treatment in the detection process are required, we have to take into account 



the parton distribution and give e s the energy dependences ||12|| . Note that all e's are 
complex numbers, i.e., e s ^™ ,d = \e s '™' d \e l 't >s ' m ' jgj. 

Then we estimate the shift of the event number induced by the new physics in ICARUS 
and OPERA, and consider the condition that the difference from the standard model be- 
comes significant. Because of low statistics and high neutrino energy, these two experiments 
are not assumed to see their energy spectra and we follow the assumption. 

We set the condition for the criterion to insist the observation of new physics that the 
deviation of the number of event caused by the new interactions, |A"^ P |, is grater the error 
for the number of event expected by standard interaction, . Here N^ M and N^ p are 

iV T SM = C J dE v f^{E v )P^ VT {E v )a UT {E v )eft{E v ), (4) 



iV r NP = C J dE u f u ^(E u ){P u ^ UT (E u ) - P^ UT (E u )}a UT (E u )eS(E u ). (5) 

Here f u is the muon neutrino flux and <j Vt is the charged current cross section for tau 
neutrino, and they are given in Ref.||. C is defined as iV pot Md e t[kton]A^A x 10 9 , where 
Na is the Avogadro's number, iVpot is the total number of proton on target in the beam 
production, and M& ct is the detector mass. Incidentally ICARUS has a mass of 5ktons and 
that of OPERA is 1.8ktons. The standard oscillation probability without exotic interactions, 
PvUut an d that with the new physics effects, P u ^ Vt can be approximated as 

P^ T = sin 2 2^cos^ 13 (^Ly, (6) 
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+ 2 sin 26 23 cos 26 23 cos 2 0i3 (i?e[e* T ] - i?e[ej r ]) 
- 2 sin 26 23 cos 2 6 13 (lm[e; T ] + Im[e* T ]) ^ 
+ 4sin 3 20 23 cos 2 ^e[ e -] (A L ) (^Jl L ) , 
where 0j,'s are the lepton mixing angles defined with the following mixing matrix 
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5m 31 is the larger mass difference, and L is the baseline length, 732km. To derive eq.(^), we 
treat Sm^ and e as the perturbations and adapt high energy approximation 5m 31 L <C E u , 
which is now appropriate assumption, and then extract terms of 0(e,8m 31 ). The detail of 



this derivation is given in Ref . [|i~2|, [T3|j . We use the detection efficiency, eS(E l> ), estimated 
in Ref.|^, The condition that we set in the beginning of this paragraph, the number of 
events induced by new physics is larger than the error in the standard oscillation assumption, 
can be described as follows: 
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There are three kinds of errors with different origins: (i) The statistical one, <r s t a , which is 
estimated by ^/ A^ M . (ii) The errors coming from the uncertainties of the oscillation param- 
eters. These are represented as the second term in the square root of eq. (|9]) according to the 
error propagation prescription, where X Q is one of n parameters included in N^ M and <J pa r a 
is its uncertainty. Exactly speaking, this treatment works well only when N^ M depends 
linearly on all the parameters. In general, dependence of A a s is not so simple. However 
since here the oscillation probability can be approximated as sin 2 2# 23 cos 4 9is(5m 2 2 L / 4E U ) 2 
in almost all energy region referred now, by regarding sin 2 2# 23 cos 4 9i 3 {5m\ 2 ) 2 as one P a ~ 
rameter included in A^ M , this method can be justified approximately. It is expected that 
the precision of sin 2 2# 23 , A(sin 2 # 23 ), becomes 1% and that of 5m 2 1; A(<5m 2 1 ), reduced to 
be 3% in the next-generation-experiments || . Moreover the uncertainty of cos 2 #13 does not 
affect the error estimation because of the smallness of #13. From these consideration the cr par 
is calculated as 

a 2 ar ~ A(sin 2 2# 23 cos 4 6 13 )(5mj 2 ) 2 

+ 2(sin 2 #23 cos 4 # 13 )(5m 2 2 ) A(5m 2 2 ). (10) 

(iii) The systematic error, a sys , which is given in Ref.[||, |5], [7|. The ICARUS collaboration 
reports the relation between the detection efficiency and the background event rate, and give 
some studies with different event selection rules in Ref. Among them, we pick out two 
cases, and refer as ICARUS-A and B. The efficiencies and errors for these cases are given in 
Table |. 





M de t eff C7 2 ys 


N^ M total error 


ICARUS-A 


5kt 0.081 11 


40.4 


7.4 


ICARUS-B 


0.047 1.5 


23.5 


5.2 


OPERA 


1.8kt 0.091 0.75 


16.3 


4.3 



TABLE I: The experimental parameters and the total error for ICARUS and OPERA. Here, we 
assume iV po t = 5 x 4.5 x 10 19 pots. Total errors are given by the left hand side of eq.([|). 

Now, we show our results. Here, we expect that the CERN proton beam will achieve 
4.5 x 10 19 pot per year and assume 5 year running. The total error in each experiment is 
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indicated in Table |. For numerical calculation we use the following theoretical parameters: 

sm0i2 = l/2, sin0 23 = l/y/2, sin0 13 = 0.1, 

5m 2 31 = 3 x 10" 3 [eV 2 ], 6m 2 2l = 5 x 10" 5 [eV 2 ], (11) 
5 = tt/2. 

Since eq.(^) is good approximation in this context, we can expect the numerical results do 
not depend on 5m 2 i, sin# 12 , 5, and sin# 13 . 

We firstly categorize the flavor-changing interactions into two classes by Lorentz and 
SU(2)l properties. The first one corresponds to the case that there exists an effective 
flavor- changing interaction of singlet/triplet type, say 

{I T T 2 Cq)(qCh 2 l^ + h.c. and/or {I T T a Cq)(qC^T a l^ + h.c. D (z/ r 7»(J 7(T w) + h.c, (12) 

where /, q are the lepton and quark doublets respectively and r a is the Pauli matrix and 
here C denotes charge conjugation. This type of an interaction is induced by the exchange 



of SU(2)l singlet and/or triplet scalar [|I2~ , T%\ . New physics effects are expected to be the 
same order of magnitude at the source, the matter and the detector for this case. In this 
case, there is a constraint for e's from SU(2)l counter part process of eq. ([l2|) like t~ — > fi~ir 
l4| . Isospin breaking effect somewhat relaxes the limit. However after all, we have to set 



the e* r < C(10" 2 ). Therefore we set |e^ r m ' d | = 0.01 for the new physics parameters. The 
left plot of Fig.|l] shows dependence of the region that satisfy the condition eq.(^). The 
inside of this contour represents that we can observe the new effect at lcr confidence level. 
The right one is a section of the left one at e™ T = O.Oli, though there the contours denote 
the difference of the event number due to the new effect from the expected number with 
standard oscillation assumption. These behaviors can be very well understood by eq. ([!]). It 
strongly depends on the phases of the exotic interaction couplings. In some regions the total 
effect can exceed the range of the error. 

Next, we assume that |efj r | = 0.01 and |ep d | = 0. This parameter set corresponds the 
situation that the new interaction is doublet type, namely, 

(l T Cd R ){qC ] n R ) + h.c. D {vr^R){dRU L ) + h.c. (13) 



This is induced by the SU (2) L doublet intermediation. From the relation |L5 



«75^ = ; <9a(M7 CT 75^), (14) 
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FIG. 1: The left plot represents the phase dependence of the contour which denotes la deviation 
of the event number from that with the standard expectation when the magnitudes of all the 
exotic couplings, |e^"' d | are assumed to be 0.01 in the case of ICARUS-A. The right one is the 
expected event number due to new interactions, A^ p in the presence of the new physics couplings 
| e iur| = 0.01 and e™ = O.Oli in ICARUS-A. This plot is almost the same as that of e™ T = 0. In the 
region around the (— vr/2, — 7r/2) and (tt/2, tt/2), the deviation from the standard case becomes 
significant beyond the error indicated in Table |. 

the doublet mediation amplitude gets the enhancement factor, m^/(m u + m^), in the pion 
decay process. On the contrary there is no such an enhancement on the propagation and 
detection processes. Therefore e* T is much bigger than ep d and hence only e^ r can contribute 
to the oscillation phenomenon. This enhancement allows us to search the smaller exotic 
coupling, which included in the elementary process, by (9(10 -2 ) than that in the singlet and 
triplet cases, eq.(|l~2"D. In this case there is a constraint on the effective coupling from the 
SU{2) L counter process, r~ — > fi~+n°, namely e* r < O(10~ 2 ). The results of this calculation 
are presented in Fig.|2|. This is also well understood by eq.([7|). In the region around <j) s = 
±7r/2 the gaps from the standard oscillation expectation become large. However, they do 
not reach la level of significant even if the events of ICARUS and OPERA are combined. 

Finally, we discuss the sensitivity for different oscillation parameters from eq. flTT|) . As we 
noted again and again, eq.(^) is a very good approximation in this context. It shows that 
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FIG. 2: Expected event number made by the new physics coupling |e 
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in ICARUS-A. The horizontal axis represents the complex phase of the e* T and the vertical axis 
indicates the event number deviation from standard model expectation, . 

the sensitivity does not depend on the sub-leading oscillation parameters, Sm^i, sin 812, 5, 
and also independent of sin^ since all terms depend only on cos^. According to the fact 
that the statistical error is dominant among the three kinds of errors, the condition eq.@ 
reduces to P v 



_ pSM > 



Pv^L, Vt - From eqs.((|) and (0), it is found that this inequality 
depends essentially on two parameters, Sm^ and sin 823. For example even if Sm^ ~ 5 x 1CT 3 
eV 2 , though the event number itself increases, we have almost same sensitivity as that with 
Sm 2 31 = 3 x 10- 3 eV 2 . Note that here Sm^L/AE < 1 |7 



To conclude, we summarize the results of the paper. We roughly evaluated the feasibility 
to search the new physics with ICARUS and OPERA, and got the following conclusions. 

• There is a possibility for the observation of the new interaction effects whose magnitude 
is |e^"' d | > O(10~ 2 ). This value is around the current bound. 

• The effects strongly depend on the phases of the couplings. 

• To explore much smaller e, more statistics is necessary. For this purpose, it is important 
to raise the efficiency even if the some background events contaminate the total tau 
events. 
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